Quantum fluctuations of the ultracold atom-molecule mixtures 
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We investigate evolution of the quantum coherence in the ultracold mixture of fermionic atoms 
and bosonic dimer molecules. Interactions are there experimentally controlled via tuning the ex- 
ternal magnetic field. Consequently, the fermionic atoms and their bosonic counterparts can be 
driven to a behavior resembling the usual BCS to BEC crossover. We analyze in some detail how 
this quantum coherence evolves with respect to time upon a smooth and abrupt sweep across the 
Feshbach resonance inducing the atom-molecule quantum fluctuations. 
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ATOM SUPERFLUIDITY 

The recent experimental techniques for trapping and 
cooling of the atomic vapors enabled exploration of the 
extremely low temperature regions where quantum ef- 
fects play a crucial role. An example can be the Bose Ein- 
stein condensation (BEC) produced out of the bosonic 
atoms in alkali metals, polarized hydrogen, etc. Phase 
transition to the BEC state is there triggered purely 
by the quantum statistical requirements which lead to 
macroscopic occupancy of the lowest energy level and 
can occur even in absence of any interactions. Recent ac- 
tivities in the field of ultracold atomic systems focus on 
application of similar techniques to fermionic atoms like 
6 Li or 40 K (besides even number of nucleons they consist 
of odd number of electrons). At ultralow temperatures 
such quantum effects like the Pauli principle play con- 
siderable role, but eventual quantum phase transitions 
would be allowed only if fermionic atoms get correlated 
via interactions. 

Interactions between trapped atoms are routinely in- 
duced by applying the magnetic field to fermion atoms 
prepared in several (two or more) hyperfine configura- 
tions. From elementary considerations 0] it turns out 
that the involved hyperfine states experience the effec- 
tive scattering described by a potential whose magnitude 
and sign depend on the applied field B. In particular, the 
various (so called) Feshbach resonances can take place. 
On this basis there was proposed a mechanism of the res- 
onance superfluidity with a transition occurring near 
the Fermi temperature T c ~ Tp, Besides the isotropic 
phase there has already been observed also the exotic 
p-wave superfluidity 0,0]. 

A unique possibility of controlling the effective inter- 
actions gives a chance for the experimental realization of 
the BCS to BEC crossover. The BCS limit corresponds 
to a case of weakly attracting fermion atoms which get 
coupled into the large Cooper pairs. In opposite limit, 
the tightly bound diatomic molecules are formed which 
ultimately can undergo transition to the BEC. Experi- 
mentalists are able to switch between these limits in a 
controllable manner and change of the interactions can 



be performed cither adiabatically (by slowly varying the 
field) [j| or in a non-adiabatic way (via the sudden sweep) 

In this short paper we investigate the quantum fluctua- 
tions induced by the time-dependent change of the inter- 
actions. We focus on a situation when the magnetic field 
is detuned from the resonant value Bq towards the far 
BCS regime at higher field B > Bq. We consider two dif- 
ferent processes: the smooth and the sudden switching. 
The fast sweep has been discussed in the literature but 
with some ambiguous conclusions concernig evolution of 
the order parameters with respect to time 0,00]. From 
our analysis we find that both parameters do oscillate in 
a damped way. 



HEISENBERG EQUATIONS 

In a close proximity to the Feshbach resonance (i.e. 
when B ~ Bq) the ultracold fermion atoms coexist and 
interact with the diatomic molecules. On a microscopic 
basis this situation can be described in terms of the two 
component boson fermion Hamiltonian Q 
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which has been known and studied in the solid state 
physics by J. Ranninger and coworkers jlOj as a phe- 
nomenological model for the high temperature supercon- 
ductivity. In the present context |QJ describes the atoms 
in two hyperfine states denoted symbolically by a =]" and 
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j. The second quantization operators c[ 
spond to fermion atoms with energy e£ = h 2 h 2 /2m and 
to diatomic molecules with energy = ?i 2 k 2 /2(2m). 
The effect of external magnetic field is included via 
the detuning parameter v which shifts the boson ener- 
gies and hence affects efficiency of the boson-fermion 
coupling g [llj. As usually /i is the common chemi- 
cal potential and we use the grand canonical ensemble 
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to ensure the conservation of the total particle number 

We are interested here in studying the time dependent 
evolution of fermion and boson occupancies together with 
the corresponding order parameters. For this purpose we 
derive the Heisenberg equations of motion which for the 
Hamiltonian are given by 
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which are identical with expressions (5) and (6) in the 
Ref. 0- One next replaces the boson operators by 
their time-dependent expectation values b(t) = (bo) and 
b*(t) = (bl). 

In the stationary case when all parameters in JHJ) are 
time independent we can derive various expressions for 
the static expectation values [ljj. Hamiltonian © has 
formally the following structure 



k • <?k + const, 
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so the pseudospin c? k behaves as though affected by 
a fictitious magnetic field = (— A', A", £k) where 
where £ k = e k - fj,, Eq = e B + 2v(B) - 2/x and we set A' + zA" = 5(60)- Following Anderson we can solve 

this problem for arbitrary temperature. In anal- 
ogy to the Weiss theory of ferromagnetism we obtain 
that a magnitude of the pseudospin expectation value 
is I (5k) I = tgh + |A| 2 /2fc s T j. Determining the 

angle between the z and xy components of the vector /i k 
we finally arrive at the stationary equations |l0| 



fi = 1 . In general equations I|2I5(I are not solvable exactly. 
In the next section we briefly discuss an approximate 
method which shall be valid for the ground state and for 
very low temperatures. 



THE SINGLE MODE APPROACH 

For temperatures close to the absolute zero we can 
neglect the excited (finite momentum) boson states. It 
is sufficient to restrict attention to the q = boson level 
because it is macroscopically occupied. In such single 
mode approach 0, the initial Hamiltonian reduces 
to 
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Following Anderson 0] we introduce the pseudospin no- 
tation 



FLUCTUATIONS OF THE ORDER 
PARAMETERS 
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such that i7 k = i (er k ± icx k ) and cr k = [cr k , cr k ] are the 
usual Pauli operators. In the single mode approach we 
rewrite the Heisenberg equations I|2I5|) using the pseu- 
dospin notation 



dt 
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In the symmetry broken state (for T < T c ) the two 
component model Q is characterized by two order pa- 
rameters: b(t) and another one of the fermion subsystem 
defined as x(T) = X)k( c - k l Ck T)- These quantities are 
complex. In the stationary case they are proportional to 
each other as can be seen from the equation l|12fl . How- 
ever, this relation is no longer valid when the Hamiltonian 
© depends on time. Evolution of the order parameters 
b(t) and x{t) with respect to time must be determined 
by solving the Heisenberg equations (|7I8|I subject to some 
boundary conditions. 
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FIG. 1: Variation of the fermion |x(t)| an< A boson \b(t)\ order 
parameters caused by the detuning from the Feshbach reso- 
nance. Initially, for t < 0, the system is in the stationary state 
with v(t<0) = fi where the magnitudes of both order param- 
eters are static. Upon lifting the molecule level there appear 
the damped oscillations. In the left h.s. panel we illustrate 
the case of an abrupt detuning v(t > 0) = 0. lg and in the right 
h.s. a smooth detuning u(t) oc t. 



We analyze here such dynamics assuming that initially, 
for t < 0, the mixed atom molecule system is exactly on 
the Fesbach resonance v = fi. We also assume that the 
boson order parameter is real and b(t < 0) = 1 so that 
the fermion order parameter is real too. Value of x(*<0) 
was determined from the equation 112|) . For simplicity 
we focus on the ground state and set the boson fermion 
coupling j as a unit or all the energies appearing in our 
study. 

For time t > we change the detuning parameter v in 
the following ways: a) via the sudden detuning as pre- 
viously discussed in the Refs 0, S and b) thr ougn 
gradually increasing v(t) — \x cx t. Avoiding any con- 
straint solutions we solved numerically the Heisenberg 
equations (|7I8|I by means of the Runge-Kutta algorithm. 

By increasing v the boson fermion system is pushed to 
the far BCS regime. For v = O.lg at t — ► oo the order 
parameters b(t) and x(t) should reach very small values. 
Figure 1 shows that indeed temporal evolution occurs 
into such asymptotics and practically is achieved after 
several oscillations. In both cases the oscillations are 
clearly damped in agreement with the previous study by 
K. Burnet and coworkers 0- However, this process of 
damping is sensitive on a particular profile of the time 
dependent detuning. This can be seen from the figure 1 
and also in the next figure 2, where we plot the phase 
9 of the boson order parameter b(t) = \b(t)\e ie ^\ For a 
smooth switching the oscillations seem not to look regular 
at all. 
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FIG. 2: Evolution of the phase 9 of the boson order parameter 
b(t) with respect to time t for: an abrupt detuning (left), and 
for a smoothly increasing detuning v(t) oc t (right). Instead 
of the bare angle 9 we plot the function sin 9 = lmb(t)/\b(t)\. 



SUMMARY 

We studied the dynamics of the ultracold fermion 
atoms upon the sudden and gradual detuning from the 
Feshbach resonance. Such situation can be experimen- 
tally realized by switching the external magnetic field 
from Bq to the higher values. From the selfconsistent 
numerical solution of the equations of motion we find 
that the order parameters start oscillating with the am- 
plitude decaying in time. Such damped oscillations de- 
pend on the specific form in which the detuning v(t) is 
carried out. We also remark that the quantum oscilla- 
tions of the order parameters turn out to be damped even 
on the level of the single mode approach without taking 
into account scattering to the finite boson momenta. 
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